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Neutrino oscillations in a hot and dense astrophysical environment such as a core-collapse 
supernova pose a challenging, seven-dimensional fiavor transport problem. To make the 
problem even more difficult (and interesting), neutrinos can experience collective oscil¬ 
lations through nonlinear refraction in the dense neutrino medium in this environment. 
Significant progress has been made in the last decade towards the understanding of 
collective neutrino oscillations in various simplified neutrino gas models with imposed 
symmetries and reduced dimensions. However, a series of recent studies seem to have 
“reset” this progress by showing that these models may not be compatible with collec¬ 
tive neutrino oscillations because the latter can break the symmetries spontaneously if 
they are not imposed. We review some of the key concepts of collective neutrino oscil¬ 
lations by using a few simple toy models. We also elucidate the breaking of spatial and 
directional symmetries in these models because of collective oscillations. 

Keywords: collective neutrino oscillations; spontaneous symmetry breaking; core-collapse 
supernova. 
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1. Introduction 

Neutrinos are abundantly produced in hot and dense astrophysical environments 
such as the early universe, core-collapse supernovae (SNe) and black-hole accretion 
discs in which they are also instrumental in the dynamical, thermal and chemical 
evolution of these environments. In a SN, for example, the hot proto-neutron star 
(PNS) at the center quickly cools down by emitting ~ 10®® neutrinos in all flavors 
within just ~ 10 seconds.^i^ Outside the PNS electron neutrinos and antineutrinos 
influence the supernova dynamics and nucleosynthesis through reactions 

i^e + np + e~, i'e + Pn + e~^. (1) 

Because the weak-interaction states (/3 = e,p,,T) and the mass eigen states 
jr'i) (i = 1, 2, 3) of the neutrino do not coincide with each other, neutrinos can mu¬ 
tate from one flavor (or weak-interaction state) into another during propagation, 
which is known as the neutrino flavor transformation or neutrino oscillations.l^The 
oscillation of neutrinos between the electron flavor and other flavors can have im¬ 
portant physical consequences in astrophysical environments through the reactions 
in Eq. Q. 
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Even in the coherent regime (i.e. without neutrino absorption, emission and 
collision) the transport of neutrino flavors poses a very challenging problem. One 
of the challenges is that, in an environment with a large neutrino flux, the flavor 
transport can become nonlinear because of the neutrino-neutrino coupling. In most 
cases the nonlinear equations of motion (e.o.m.) which govern neutrino oscillations 
have to be solved numerically. The computational investigations by various groups 
show that a dense neutrino medium can experience collective oscillations during 
which neutrinos of different momenta oscillate cooperatively.^^®^ To obtain theo¬ 
retical insights into these numerical results, several toy models such as the bipolar 
modeP^®^ have been proposed which can be solved analytically and which bear 
some of the qualitative features of the more realistic models. Linear stability anal¬ 
ysis provides another useful tool to predict the physical regimes where collective 
neutrino oscillations may or may not occur .1^ 

Another challenge in solving the problem of neutrino transport is simply due 
to its large dimensionality: the neutrino transport in any real astrophysical envi¬ 
ronment even without flavor oscillations involves 1 temporal dimension, 3 spatial 
dimensions and 3 momentum dimensions. This challenge has been largely bypassed 
in the previous studies which are limited to the models with imposed symmetries. 
For example, the isotropic and homogeneous condition is usually assumed for the 
early universe. However, even though such symmetries can exist in the e.o.m. and 
may also exist in the physical systems at some point, they may not be preserved 
during neutrino oscillations.^^ Indeed, a series of recent work suggest that these 
symmetries can be broken spontaneously by neutrino oscillations 

The goals of this paper are to review some of the key concepts of collective neu¬ 
trino oscillations through the bipolar neutrino gas model and to explain why spatial 
and directional symmetries can be broken spontaneously by collective oscillations. 
We will leave out many important topics such as three-flavor oscillation£3111[2ll 
some of which can be found in an earlier review.®! 

Throughout this paper we adopt the natural units with h, = c = 1. We also 
assume that neutrinos are relativistic with speed u ^ 1. 

2. Collective neutrino oscillations 
2.1. Equation of motion 

Assuming the validity of the mean-field theoryjSlEll we will use the (Wigner- 
transformed) neutrino (flavor) density matrice^^ Pp(t, x) and Pp(t,x) to describe 
the flavor content of the neutrino and antineutrino of momentum p at time t and 
position X. For simplicity we will consider neutrino oscillations between the elec¬ 
tron flavor and another active flavor, say r. In the weak-interaction basis (or flavor 
basis). 


■ Pt 
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■ n Pv' 

.(Pp")* pV. 
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.(Pp")* pV. 
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The diagonal elements of the density matrices in the above equation are propor¬ 
tional to the number densities of the neutrino or antineutrino in the corresponding 
flavors and momentum state, and 


mg 




(3) 


are the total number densities of the neutrinos and antineutrinos in flavor j3 {/3 = 
e,T), respectively. The off-diagonal elements of the density matrices contain the 
information of flavor mixing. In the rest of the paper we will use the weak-interaction 
basis exclusively unless otherwise stated. 

In the coherent regime a neutrino can change its flavor through refraction when 
it passes through a dense medium. In the lowest-order and without gravitational 
redshift, the e.o.m. that govern the neutrino flavor transformation ar^^^HaSl 


T V • V)pp — i[Hvac T T Pp], (4^) 

(l9t -|- V • V)pp = —i[—Hvac + Hmat + Pp]) (4b) 


where v = p/E is the propagation velocity of the neutrino with E 
energy of the neutrino. The vacuum Hamiltonian is 

_ Am^ — cos 20V sin 29^ 

4i? sin 20V cos 20V ’ 


IpI being the 


(5) 


where Am^ is the neutrino mass-squared difference, and 0v is the vacuum mixing 
angle within range (0,7r/4]. Here Am^ > 0 and Am? < 0 correspond to the nor¬ 
mal neutrino mass hierarchy (NH) and the inverted hierarchy (IH), respectively. 
The matter potential in Eq. Q arises from the coherent forward scattering of the 
neutrino by the charged leptons in the medium through the charged-current weak 
interaction.!^ When the densities of p, and r leptons are negligible, the matter 
potential can be written as 


Hmat = V2 Gf 


Ue 0 
0 0 


^(1 + f ^ 3 ), 


( 6 ) 


where Gp is the Fermi constant, Ue is the net electron number density, A = v^GFU-e, 

I is the identity matrix, and is the third Pauli matrix. The neutrino(-neutrino 
coupling) potential in Eq. Q stems from the coherent forward scattering of the neu¬ 
trino in question by ambient neutrinos through the neutral-current interaction,!^^!!^ 
and it can be written as 

= V2Gf J ^^^(1 - V • V )(0p' - PpO) (7) 

where the physical quantities with primes are for the ambient neutrinos. 

We will be interested in the regimes with large matter densities where the (ef¬ 
fective) neutrino mixing angle in matter is small and 


%/2GFne»|Am^/2£;| 




(A-po;)^, 


— Hyac H" Hjnat 


( 8 ) 
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where t} = +1 and —1 for NH and IH, respectively, and 

a;(i?) = ^^cos20v. (9) 

Zh/ 

We have ignored the trace terms here which have no impact on neutrino oscillations. 
For antineutrino, 

Btj = —Hvac + Hmat - > (A + TJUl) — = (10) 

Therefore, for the purpose of neutrino oscillations one can treat an antineutrino of 
oscillation frequency w and energy E as a neutrino with oscillation frequency —lo 
and energy —E. We will adopt this convention throughout this paper. 

In the absence of neutrino absorption, emission and collision, it is useful to 
define (reduced) neutrino flavor (density) matrices 

Puj>0,v CC Pp, Pu;<0,v C)C Pp (f f) 

with normalization condition 

trpcd.v = 1- (12) 

The diagonal elements of Poj,v give the probabilities for the corresponding neutrino 
(uj > 0) or antineutrino (w < 0) in the corresponding weak-interaction states. 
Neutrino flavor matrices obey the e.o.m. 

(a* + V • V)p^^-o = -i[B,^ + (13) 


2.2. Bipolar model 

The study of collective neutrino oscillations begins with the bipolar model.^l^lfl^This 
model describes a homogeneous and isotropic system with Vg and populations of 
single energy Eq at time t — 0. The e.o.m. of neutrino oscillations is 


idtPu: 


(A-r?w)y + H, 


Qlo 


with CO = iwo, and the neutrino potential is 


(14) 


— p{PbjQ 0:Q—ujQ')t 

where 

Tip 

a = — 

Ui, 

is the ratio of the density of the antineutrino to that of the neutrino, and 

pL = V^Gprii, 


(15) 

(16) 

(17) 


is a measure of the strength of the neutrino-neutrino coupling. 


If the neutrino density is negligible, Eq. (14) has solution 

e±u,o{t) ~ exp |^-i(A =F »?wo)tyJ £'±t.uo(0) exp |^i(A =F ?7‘^o)iy 


(18) 
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Because 

^?iha;o(0) 

are diagonal and commute with one has 


1 0 
0 0 


1 0 
0 0 


(19) 


( 20 ) 


In other words, neutrino oscillations are suppressed in the presence of dense matter 
if there is no ambient neutrino. 

However, the presence of dense matter does not necessarily suppress collec¬ 
tive oscillations in the dense neutrino medium. This can be seen by performing a 
“corotating-frame” transformation to the neutrino flavor matrixl^^ 


~Qu.it) = 


( 21 ) 


which removes the time evolution due to matter. Under this transformation Eq. (14) 
becomes 

0-3 


idt gu.it) = 


-rfuj- 


H,. 


Quj 


( 22 ) 


where = fJ-igujo~(^Q-ujo)- Eqs. (14) and (22) are equivalent, and p and phavethe 
same diagonal elements. Therefore, if there exists collective flavor transformation 
in the bipolar model, a uniform dense matter will not suppress itQln the rest of 
the paper we will always work in the corotating frame by taking A = 0, but we will 
drop the tilde for simplicity. 

To see that a bipolar system can experience collective flavor transformation we 
note that Eq. (22) is invariant under global phase transformatiorl^ 


gu.it) 


3-i0CT3/2 


gu.it) 


, 100 - 3/2 


(23) 


where ^ is a constant. This symmetry leads to the conservation of the “flavor lepton 
number (density)”!^ 


C = n^tT[igu.o - ag-u.o)a-3] = (n^^ - - (rip^ - np^). 


(24) 


The phase symmetry in Eq. (231 also implies the existence of a collective precession 
solution: 


g^{t) = 


(25) 


where H is a common oscillation frequency independent of the oscillation frequency 
uj of the individual neutrino. 

The bipolar model can be solved analytically.^ It has also been shown that 
the bipolar model is equivalent to a gyroscopic pendulum in flavor space,^ and the 
collective precession solution in Eq. (251 corresponds to the precession motion of this 


^However, collective neutrino oscillations can be suppressed by a very large matter density in SNe 
as shown in Ref. 1401 
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pendulum. However, this collective precession solution can not be achieved exactly 
for arbitrary initial conditions just as a gyroscopic pendulum usually experiences 
both nutation and precession. 

The bipolar model is, of course, very different from the neutrino medium in a 
realistic astrophysical environment. However, it does provide explanations for some 
of the important results of collective neutrino oscillations in more sophisticated 
models. For example, one can show that in the collective precession solution neu¬ 
trinos will remain in the weak-interaction states if the neutrino density is larger 
than a critical valu^^^ [see the discussion in Sec. 2.4 and around Eq. ©]■ This 
corresponds to the “sleeping-top” regime of a gyroscopic pendulum where the fast 
spinning top or pendulum defies gravity and spins without wobbling. When the 
neutrino number density decreases below the critical value, the flavor pendulum 
begins to wobble about the precession solution. Similar phenomenon was indeed 
observed in numerical simulations of neutrino oscillations in SNe (see Fig. [^. 




Fig. 1. Left panel: The collective precession solution of a bipolar model with time-varying neutrino 
density nu (thick lines) and the numerical solution for this model (thin lines). The dashed and 
solid lines represent cosi?i = tr(pQjQ(T 3 ) and cos»92 = —tr(p—i^o^'a), respectively. Right panel: 
The energy-averaged neutrino survival probabilities for three representative neutrino trajectories 
as functions of radius r in a spherical supernova model. The figures are adapted from Fig. 3 of 
Ref. |43| and Fig. 2 of Ref. respectively. 


2.3. Spectral swap/split 


The collective precession solution in Eq. (251 also exists in a homogeneous and 
isotropic neutrino gas with continuous neutrino energy spectra.^ We will again 
assume that the neutrinos and antineutrinos are in the weak-interaction states at 
time t = 0. Because Eq. ([^ is invariant under p —> —p (without changing we 

can use 


1 0 
0 0 


Quit = 0 ) 


(26) 
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as the initial condition for all neutrinos and antineutrinos by replacing the neutrino 
energy spectra with a single effective spectrum 


g{u}) cx 


dE 


dw 


r tr(p_EO-3)t=o if w > 0, 
\ -tT{pE<J3)t=0 if w < 0 


with normalizatior0 


poo 

/ g(w)da; = l. 
Jo 


(27) 


(28) 


The e.o.m. (22) also applies to the neutrino gas with continuous neutrino energy 

(29) 


spectra except that the neutrino potential is now 

/ OO 

g{uj)Qu, dw, 

-OO 

where 

g = V2GF{n^^ - n^^)t=o- 


The definition of a in Eq. (16) can be generalized to 

/°oo 9 {^) 


(30) 


(31) 


Obviously, the bipolar model is a special case of the continuous-spectrum model 
with 


g{ui) = S{uj — ujq) — aS(uj + ojo)- 


(32) 


If the neutrino gas is in the collective precession mode initially and expands adi- 
abatically afterwards, it should remain in the collective mode. However, according 
to Eq. (22), when the neutrino density becomes negligible (at time t > ti), neutri¬ 
nos of different energies should oscillate independently with their own oscillation 
frequencies: 




(33) 


This apparent contradiction is resolved if guj{ti) is diagonal. Being diagonal does 
not necessarily mean that Qui(ti) = £itj(0) = diag[l, 0]. For example, if the neutrino 
gas has a positive flavor lepton number C [see Eq. ( p4| )] and a bipolar-like spectrum 
[i.e. g{uj) is positive (negative) if w > 0 (w < 0)] and if the neutrino mass hierarchy 
is inverted, the final spectrum can be split at the critical energy 


Eq 


I Am^ 

2Hn 


(34) 


Other equivalent normalization conditions also exist in the literature. The normalization con¬ 
dition in Eq. l|28|l is used in anticipation that the Ue flux is the largest among all flavors in a 


SN. 
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such that 


rdiag[l,0] if w > r^o, 
\diag[0,l] if w < rio, 


(35) 


where fig is the collective oscillation frequency at /x = cP and can be determined 
by the conservation of The above result implies that Vf, and Vt have swapped 
their energy spectra at energies E > Eq. This spectral swap/split phenomenon is 
commonly observed in the numerical simulations of neutrino oscillations in SNe (see 
Fig.[|). If the neutrino gas has a more complicated spectra with multiple spectral 
crossings where g{uj) = 0, the final spectra can exhibit multiple spectral splitsl^^ 



Fig. 2. Left panel: The split neutrino spectrum (dotted line: fully adiabatic; thick solid line: 
numerical) produced from a box-like initial spectrum (thin solid line) in a homogeneous and 
isotropic neutrino gas by collective neutrino oscillations, where u! = |Am^|/2i?y and = 

tr(gi,j). Right panel: The final neutrino survival probability as a function of neutrino energy Ei, 
and neutrino emission angle i?o ou th® surface of the proto-neutron star in a numerical calculation 
employing a spherical supernova model. The figures are adapted from Fig. 1 of Ref. |42| and Fig. 3 
of Ref. respectively. 


2.4. Flavor instability 

Although the bipolar model can be solved analytically and provide some qualitative 
understandings of collective oscillations in neutrino media, numerical computations 
are usually required to obtain quantitative results of neutrino oscillations in more 
sophisticated models. However, if the neutrinos and antineutrinos are initially in 
the weak-interaction states, linear stability analysis can be utilized to predict the 
physical regimes where collective oscillations can occur without performing large- 
scale numerical simulations. This is especially useful if a large parameter space 
needs to be surveyed because of the uncertainty in, e.g., the neutrino fluxes in 
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SN simulations. The linear stability analysis method was explained thoroughly in 
Ref. 20 Here we will use the bipolar model to demonstrate the essence of this 
method. 

We assume that no significant flavor oscillations have occurred at time t such 
that 


QbJo 


1 e 
e* 0 


e-u 


1 e 
e* 0 


(36) 


where |e| ~ |e| <C 1. Keeping only the terms up to 0{e) we can rewrite Eq. (22) in 
terms of e and e as 


idt 


e 

e 


—r]UJ — fxa 

afi 


e 

= A- 

e 


TfUJ + fl 


e 


e 


(37) 


Note that, in the absence of the neutrino medium (i.e. /r = 0), e and e are decoupled 
and oscillate with frequencies —tjuj and t]uj, respectively. In the neutrino medium, 
however, e and e are coupled, and each eigenvalue H of matrix A corresponds to a 
collective oscillation mode: 


eit) « e(t) « (38) 

where [Q, QV is the corresponding eigenvector of A. 

If H is real, e and e will oscillate and their amplitudes will remain small. If H is 
complex with 


K = Im(H) > 0, 


(39) 


the bipolar model has a flavor instability: the amplitudes of e and e will grow 
exponentially which can lead to (large-magnitude) collective oscillations. Because 
A is a real matrix, when H is complex, Vi* is also an eigenvalue of A. In this case, 
the solution with positive k will dominate over time. Therefore, for the purpose of 
flavor stability analysis, it is sufficient to find out the regimes where A has complex 
eigenvalues. 

It is straightforward to show that A has complex eigenvalues only when 

(I - a)^/r^ -I- 4cLip -I-4(1-I- a)/r?7Wo < 0 (40) 


or 


2wo 


< -w < 


2uJo 


(41) 


(1 + V^)2 ^ "^'( 1 - V ^) 2 - 

The bipolar model (with i/g and Pe initially) can have flavor instability only in IH. 
The upper limit of /r where collective oscillations can occur indeed agrees with the 
result plotted in Fig. [^a). It is interesting to note that, if the bipolar model has 
and Pr initially instead of i/e and Pe, it will have flavor instability only in NH 
because ^ < 0 [see Eq. (30)]. 

Similarly, for the continuous-spectrum model, we assume 


Qu) 


I 


0 


(42) 
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the e.o.m. (221 can be written as 

/ OO 

g{io')eu,' dw'. 

-OO 

Using the ansatz 

one obtains 

1 

[f2 — (1 — a)fi + r]Uj]Qaj = const. or Qi^ cx 


U — (1 — a)fi + rjUJ 


(43) 

(44) 

(45) 


In order for the above solution be consistent with Eq. (43), U must satisfy the 
consistency condition 

_ 1 ^ gjuj) 

M J_^n- {l-a)g. + guj 


dw. 


(46) 


Because U = U — (1 — a)g. and U have the same imaginary part, it is sufficient to 
solve 

1 r 5(^) 




/—OO H“ TjUJ 


duj 


(47) 


for U for the purpose of flavor stability analysis. A complex implies the existence 
of a flavor instability. 

Flavor stability analysis can also be applied to more complicated models. Al¬ 
though a powerful tool, the linear stability analysis method depends on the validity 


of assumption (42) and is limited to the linear regime. When there is a flavor insta¬ 


bility, numerical simulations are still needed to follow neutrino oscillations in the 
nonlinear regime. 

3. Spontaneous symmetry breaking 


The e.o.m. (13) involves seven dimensions and has not been solved in its complete 


form. Instead, various symmetries have been assumed to reduce the dimensionality 
of the problem so that a numerical or analytically solution can be found. Three 
classes of symmetries are commonly employed: the time translation symmetry that 
makes the problem time-independent, spatial symmetries that reduce the spatial 
dimensions, and directional symmetries that reduce the momentum dimensions. For 
example, a commonly used model for the early universe assumes complete (spatial) 
homogeneity and (directional) isotropyEEH These assumptions reduce the neutrino 
transport in the early universe to a two-dimensional problem. Another example 
is the SN neutrino Bulb modeP which assumes the time translation symmetry, 
the (spatial) spherical symmetry about the center of the SN, and the (directional) 
axial symmetry about the radial direction. These assumptions reduce the neutrino 
transport in the SN to a three-dimension problem. 
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However, the full e.o.m. (13) allows solutions that may or may not have these 
symmetries. If a symmetry-breaking solution becomes unstable, collective oscilla¬ 
tions can break the corresponding symmetry spontaneously even if such an (ap¬ 
proximate) symmetry exists in the neutrino medium initially.^ We will first use 
the two-beam model to illustrate the basic idea of spontaneous symmetry break¬ 
ing by collective oscillations. We will then look at the breaking of directional and 
spatial symmetries separately. 


3.1. Two-beam model 

We consider a model in which all neutrinos are emitted from the x axis and prop¬ 
agate only in the x-z plane. We assume that neutrinos are emitted in only two 
directions: 


vc = K,0,u^] iC = L,R), 


(48) 


where 0 < < 1 and ua = —up = \/l — v\. We also assume that every point 

on the X axis emits in each direction neutrinos and antineutrinos of single energy 
i?o with intensities jv and jv = o-jv, respectively. We further impose the time 
translation symmetry and the translation symmetry along the x axis. With these 
conditions Eq. (131 reduces to 

iV(; • = 


cr3 


In the above equation the neutrino potential is 


/ OO 

g{w)Q^'X' 

-OO 




aioX' 


(49) 

(50a) 

(50b) 


where 


fj. = ^2(1 - Vi • VR)GFjn, 


(51) 


and C = R,L are the opposites of C- 

The two beam model has been solved analytically.^^ Here we use the linear 
stability analysis method to demonstrate that the left-right symmetry {L <->■ R) can 
be broken spontaneously.^ We again assume that the neutrinos and antineutrinos 
are in the (almost pure) electron flavor: 


'1 ec 


'1 ec' 

.^c 

1 ~ 

ej 0 


We then define 

± El ± e_R 


El ± tR 


(52) 


( 53 ) 
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Keeping only the terms up to 0{e) we can rewrite Eq. (49) as 


\d. 


■e+‘ 


■e+' 



■a+ 



e~ 


A_ 


e~ 

_e~ _ 


_e~ _ 


where 




A- = V 


-1 


—r]U}o — cefj, a/j, 

rjujo + fJ,_ 

-rjOjQ + (2 - a)fi -a^ 

fi ■qojQ + (1 — 2a)fi_ 


(54) 


(55a) 

(55b) 


and the minus modes 


evolve indepen- 


Therefore, the plus modes 
dently in the linear regime. 

The two-beam model is a 3-dimensional model with 1 spatial dimension (z) and 
2 momentum dimensions (w, 0- Its e.o.m. (49) is invariant under the simultaneous 
interchange of the left and right neutrino beams. When this left-right symmetry 
is artificially imposed, Eq. (49) reduces to the e.o.m. of the 2-dimensional bipolar 
model [Eq. 


]. The bipolar model permits only the symmetry-preserving solutions 
which can be unstable only in IH (see Sec. |2.4[). However, because 


-I-2(1 - a)/rl, (56) 

the minus modes in the two-beam model are unstable in the same range of ^ in NH 
as the plus modes in IH. In the regime where the minus modes become unstable, 
the left-right symmetry in the two-beam model is broken spontaneously, and the 
neutrino oscillations in the two-beam model and the bipolar model are qualitatively 
different. 


3.2. Breaking of the directional symmetry 


The breaking of the left-right symmetry in the two-beam model is an example of 
how the directional symmetry can be broken during neutrino oscillations. The origin 
of this symmetry breaking lies in the current-current nature of the weak coupling. 
This is best illustrated using the bipolar model in Sec. |2.2| but without imposing 
the isotropic condition.!^ The e.o.m. of this model is 


where 




^^3 , u 


(57) 


H 


h'l' 


72 

47r 


GpUi/ 


^bv/(l V • V '){Qujq^-v' ^Q—UJq,v'') 


(58) 


with dr<;./ being the differential solid angle around direction v'. We will again as¬ 
sume that 


QuJO,^ 


1 

.4 0 . 


Q—U)q ,V 


1 Cir 

.4 0. ■ 


(59) 
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We now switch to the basis of spherical harmonics y/,m(v) and define 

ei,m = j dhv, = J dTi,. 


Using identity 


1 — V • V = 47r 


^0,o(v)Fo:o(v ) - J E 


m=0,±l 


one can rewrite Eq. (571 in the linear regime and in the spherical basis as 


where 


with 


idf 


^l = {1- a)no\ + 


^l,m 

« A/ • 

^l,m 





-rjujQ - fj.1 a^i 
-fj.1 r]ujo + afj,i 


( v^Gpn^ if Z = 0, 

fii = I -{V2/3)GFn, a 1 = 1, 

[ 0 otherwise. 


(60) 


(61) 


(62) 


(63) 


(64) 


We note that the coupling coefficients /c; of the monopole mode {I = 0) and the 
dipole modes {I = 1) have opposite signs. This sign difference is derived from the 
current-current nature of the weak coupling which is proportional to (1 — v • v^) 


[see Eq. (61)]. 

The homogeneous neutrino gas model is 4 dimensional with 1 temporal dimen¬ 


sions and 3 momentum dimensions. Its e.o.m. (57) is invariant under a simultaneous 


but arbitrary rotation of the momenta of all the neutrinos and antineutrinos. When 


this rotation symmetry in momentum space is artificially imposed, Eq. (57) reduces 


to the e.o.m. of the 2-dimensional bipolar model. The bipolar model permits only 
the symmetric solutions which have flavor instabilities only in IH. However, because 
< 0, the dipole modes in the homogeneous gas model have flavor instabilities 
in NH [see the discussion around Eq. ( [4l| )] which is clearly shown in Fig.[^ In the 
regime where the dipole modes become unstable, the rotation symmetry in the ho¬ 
mogeneous gas model is broken spontaneously which makes it qualitatively different 
from the bipolar model. 

The above result can also be extended to SN models. As mentioned previously, 
a commonly used SN model in collective neutrino oscillations is the Bulb model 
which has the (spatial) spherical symmetry about the center of the SN and the 
(directional) axial symmetry about the radial direction. We define moment basis 
function 


$m(</5) = 




(m = 0,±1,...), 


(65) 
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T T 

Fig. 3. The exponential growth of the amplitude of the monopole mode (|qo,o| esc |eo,o|) in the 
inverted neutrino mass hierarchy (left) and the dipole modes (|qi,o| oc lei^oN |qi,si| ^ ~ 

ei,_i)/2i| and |qi,ci| cc |(ei,i + ei,_i)/2|) in the normal neutrino mass hierarchy (right) in a 
homogeneous neutrino gas model. The figure is adapted from Fig. 1 of Ref. |24[ 


where (p is the azimuthal angle of v about the radial direction. Using identity 


1 — V • v' = 27r 


(1 - COSl9cOSl9')®o(‘/5)®o(‘^0 - 1 Sini?sim9' 


m—dzl 


( 66 ) 

one sees that the coupling coefficients of the |m| = 0 and 1 modes have opposite 
signs, where'd is the angle that v makes with the radial direction. Similar to the 
homogeneous neutrino gas mdoel, only the m = 0 modes can exist in the Bulb model 
because of its imposed axial symmetry. As a result, collective neutrino oscillations 
occur only in IH if the neutrino fluxes are bipolar-like. If the axial symmetry is not 
imposed, however, collective oscillations can also occur in NH through the collective 
modes with m = ±1 which break the axial symmetry spontaneously. (See Ref. 21 


for a detailed analysis and Ref. 22 for the numerical confirmation of this analysis.) 


3.3. Breaking of the spatial symmetry 

Spatial symmetries can also be broken by collective neutrino oscillations. This can 
be illustrated by using the (neutrino) Line modeP^ which is a generalization of the 
two-beam model in Sec. |3.1[ Compared to the two-beam model, the Line model 
does not possess the translation symmetry along the x axis. But for computational 
purpose we will impose a periodic condition 

Qu>x{xjZ) = q^^q{x + L,z), (67) 

where L is constant. We will again assume that neutrinos and antineutrinos are in 


the (almost pure) electron flavor [see Eqs. (52) and (53)]. It is convenient to work 
in the Fourier basis and define 




I 


e-i™'=“^e=^(a;,z) dx, 




{z) = 


e-*'"'=«"^e=^(x,z) dx, 


( 68 ) 
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where feg = 2tt/L. In the linear regime the Fourier modes with different m values 
are decoupled, and the linearized e.o.m. are 


where 



-1 


—lywo — a/i CK/i Tftq 

—fi rjojo + n 0 

mq 0 —rjcoo + (2 

0 mq fj, 


0 

mq 

a)kt —a/i 

rjuJo + (1 — ‘2a)qi_ 


(69) 


(70) 


with q = (21^/L')^J\ — v^. 

The Line model is a 4-dimensional model with 2 spatial dimensions {x, z) and 2 
momentum dimensions {uj, Q). Its e.o.m. is invariant under the translation symmetry 
along the x direction. When this symmetry is artificially imposed, the Line model 
reduces to the 3-dimensional two-beam model where only the homogeneous modes 
(with TO = 0) are allowed. The Line model allows both the homogeneous and 
inhomogeneous (to ^ 0) modes, and the latter break the translation symmetry 
spontaneously when they become unstable. 

For given values of \x and a one can find out if there is a flavor instability by 
computing the eigenvalues of A. It turns out that the (two-beam) Line model has 
the same flavor unstable regions [in terms of (/r, to)] for both IH and NH. Fig. 
shows the flavor unstable region in the (/r, to) space for the case with a = 0.8 and 
L = 207r/a;o. Note that the inhomogeneous modes with large |to| values occur on 
small distance scales, and they can be unstable at larger neutrino densities than 
the homogeneous mode does. 


4. Summary and outlook 

Neutrino flavor transport in hot and dense astrophysical environments is a difficult 
seven-dimensional nonlinear problem. Using the simplest neutrino gas model, i.e. 
the bipolar model, we reviewed some of the key concepts of collective neutrino 
oscillations. By applying the method of linear stability analysis we illustrated that 
how the directional and spatial symmetries may be broken spontaneously during 
collective neutrino oscillations. 

The realization of the spontaneous symmetry breaking commences a new stage 
of the study collective neutrino oscillations. All the previous research is limited 
to the models with imposed symmetries which have smaller dimensions than real 
physical systems do. The studies of these simplified models have yielded important 
insights into the intriguing phenomena of collective neutrino oscillations, and some 
of the toy models shall still prove useful in understanding the results in more realistic 
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a=0.8 



IjI(oo 


Fig. 4. The flavor instabilities of the Line model with a = 0.8. The color scale represents 
the largest exponential growth rate of the collective oscillation modes. The size of the periodic box 
along the x-axis is taken to be L = 20'7r/a)o, and the propagation directions of the neutrinos are 
given by unit vectors \va:,Vz\ = [ib\/3/2,1/2]. The results are independent of the neutrino mass 
hierarchy. The figure is adapted from Fig. 1 of Ref. |26| 


models. However, the incompatibility of the reduced-dimension models and collec¬ 
tive neutrino oscillations clearly points out that these models are insufficient and 
that dimensionality really matters in collective oscillations. For example, a previ¬ 
ous comparison between the results obtained using the single-angle and multi-angle 
SN models shows that neutrino oscillations in these two models of different dimen¬ 
sions can have completely different yields of heavy elements in SN neutrino-driven 
wind.l^ The results of the Line model suggest that neutrino oscillations can occur 
at larger neutrino densities in the models of multiple spatial dimensions than in the 
corresponding models with only one spatial dimension. Although a very challenging 
task, new and realistic multi-dimensional neutrino gas models must be constructed 
and employed before we can know with any confidence what true physical impact 
neutrino oscillations may have on astrophysical environments with dense neutrino 
media such SNe and black-hole accretion discs. 
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